Mon. Not. R. Astron. Soc. 000, 1-?? (1994) Printed 2 February 2008 (MN LJOgX style file vl.4) 



O 
O 



Origin of the metallicity distribution of the NGC 5128 
stellar halo 



K. Bekki/ William E. Harris^'^, and Gretchen L. H. Harris^'^ 

^School of Physics, University of New South Wales, Sydney 2052, NSW, Australia 
^Department of Physics & Astronomy, McMaster University, Hamilton ON L8S ^M_Z, Canada 
^ Department of Physics, University of Waterloo, Waterloo ON N2L 3G1, Canada 

* Visiting Fellow, Research School of Astronomy & Astrophysics, Australian National University, Weston, ACT 2611, Australia 



Accepted Received in original form 2001 



in 

(N 



ABSTRACT 



o 

Q 

■ Recent Hubble Space Telescope photometry in the nearby eUiptical galaxy NGC 

psj ' 5128 shows that its halo field star population is dominated by moderately metal-rich 

stars, with a peak at [m/H] ~ —0.4 and with a very small fraction of metal-poor ([m/H] 
< —1.0) stars. In order to investigate the physical processes which may have produced 
I this metallicity distribution function (MDF), we consider a model in which NGC 

in , 5128 is formed by merging of two major spiral galaxies. We find that the halo of an 

' elliptical formed this way is predominantly populated by moderately metal-rich stars 

with [m/H] —0.4 which were initially within the outer parts of the two merging discs 
and were tidally stripped during the merger. To match the NGC 5128 data, we find 
' that the progenitor spiral discs must have rather steep metallicity gradients similar to 

, the one defined by the Milky Way open clusters, as well as sparse metal-poor haloes 

(5% or less of the disc mass). Very few stars from the central bulges of the spiral 
galaxies end up in the halo, so the results are not sensitive to the relative sizes (bulge- 
to-disc ratios) or metallicities of the initial bulges. Finally, we discuss the effects on the 
Q ■ globular cluster system (CCS). The emergent elliptical will end up with metal-poor 

| H \ halo clusters from the original spiral haloes, but with moderately metal-rich halo stars 

C/3 i from the progenitor discs, thus creating a mean offset between the MDFs of the halo 

' stars and the GCS. Remaining questions yet to be answered concern the total size of 

the GCS population (the "specific frequency problem") and the observed existence of 
, , metal- rich globular clusters in large numbers in the outer haloes of giant ellipticals. We 

' also discuss possible differences in the MDFs of stellar haloes of galaxies of different 

Hubble type. 



Key words: galaxies: abundances — galaxiesrelliptical and lenticular, cD — galax- 
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1 INTRODUCTION 

Physical properties of metal- deficient stellar haloes in galax- 
ies are considered to provide vital clues to the understand- 
ing of early dynamical and chemical evolution of galaxies. 
In particular, the detailed investigation of structural, kine- 
matic, and chemical properties of such a "fossil record" com- 
ponent (i.e., stellar halo) of the Galaxy has revealed a possi- 
ble scenario as to how the Galaxy has developed its dynami- 
cal structures such as bulge, thin and thick discs (e.g.. Free- 
man 1987; Majewski 1993; van den Bergh 1996). Because of 
the observational difficulties in revealing three dimensional 
structure and kinematics of halo stars (with respect to their 



host galaxy) in each of the Local Group galaxies other than 
the Milky Way, the metallicity distribution function (here- 
after MDF) of halo stars in these galaxies has served to 
give some constraints on the past star formation and chem- 
ical evolution histories of these systems (Mould & Kristian 
1986; DurreU, Harris, & Pritchet 1994, 2001; Pritchet & 
van den Bergh 1998; Christian & Heasley 1986: Reitzel, 
Guhathakurta, & Gould 1998; Grillmair et al. 1996; Han 
et al. 1997; Martinez-Delagado & Aparicio 1998). The stel- 
lar halo of M31 appears to be dominated by a moderately 
high-metallicity population with [m/H] ~ —0.5 (cf. the ref- 
erences cited above), a strikingly different enrichment level 
than in the Milky Way stellar halo, and a strong indicator of 
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a different formation history (e.g. Durrell, Harris, & Pritchet 
2001). 

Although there have been many MDF studies of halo 
stars in Local Group members, galaxies beyond the Lo- 
cal Group have only recently been investigated through 
J/S'r/WFPC2 photometry and for only a small number of 
cases: the giant E/SO NGC 5128 (Soria et al. 1996; Harris, 
Harris, & Poole 1999 hereafter referred to as HHP; Harris 
& Harris 2000, HHOO; Harris & Harris 2002, HH02; Mar- 
leau et al. 2000), the edge-on SO NGC 3115 (Elson 1997; 
Kundu & Whitmore 1998), and two dwaxf ellipticals in the 
M81 group (Caldwell et al. 1998). FYirthermore, Rejkuba 
et al. (2002) have demonstrated the existence of a signifi- 
cant intermediate-age AGB population in parts of the NGC 
5128 halo, based on deep U , V, and Kg color-magnitude and 
color-color diagrams of halo stars resolved by VLT. 

A remarkable result emerging from the HST/WFPC2 
photometry of the old-halo red giant stars in NGC 5128 is 
the dominance of moderately metal-rich stars in the range 
-1 < [m/H] < 0.0 (HHP, HHOO, and HH02). HHOO and 
HH02 suggested that the relatively high mean abundance 
and small fraction (~ 10%) of metal-poor stars with [m/H] 
< — 1 can give strong constraints on the total mass of dwarf 
galaxies that could be accreted by NGC 5128 and then dis- 
rupted to form its faint stellar halo. 

To date, however, there are only a few attempts at phys- 
ical interpretation of this material. HHOO discussed a one- 
zone model of chemical evolution of this galaxy and sug- 
gested that it experienced two fairly distinct stages of halo 
formation: an early "accreting box" stage when the galaxy 
is assembling through infall of star-forming gas clumps, and 
then a major stage of "closed-box" evolution when the star 
formation proceeds in the galaxy with little gas infall. HH02 
extended this basic picture further to a more generalized 
accreting-box formation model, in which the infall rate of 
unenriched gas is envisaged to start at a high rate and 
then dies away exponentially, while star formation continues 
throughout until the gas supply is exhausted. With appro- 
priate choices of gas infall rate and the effective chemical 
yield, excellent overall matches to the observed MDFs can 
be obtained. Recently, Beasley et al. (2002a, b) have fur- 
ther discussed the origin of the stellar halo MDF within the 
context of a semi-analytic model of galaxy formation (GAL- 
FORM), based on the Cold Dark Matter (CDM) picture. In 
this model, star formation can take place in two modes, a 
"quiescent" mode within the relatively unenriched pregalac- 
tic clouds, and a "starburst" mode whenever two roughly 
equal clouds suddenly merge. They demonstrated that the 
global metallicity distribution within model galaxies com- 
parable to NGC 5128 in size is broadly consistent with the 
observed MDF, although the model is unable to provide any 
information on the spatial distribution or radial change in 
the MDF. 

These previous studies are based on one-zone models 
of chemical evolution and thus cannot spatially resolve the 
various distinct stellar components (i.e., halo, bulge, and 
disc) in a galaxy. Accordingly, comparisons of these models 
with the MDFs observed at different places in NGC 5128 
have limited value. Numerical simulations which enable us 
to resolve spatially the halo, bulge, and disc components 
and thereby to investigate the MDF for each of these com- 



ponents should be very helpful for better understanding the 
formation of galactic stellar haloes. 

The purpose of this paper is to investigate the origin of 
the MDF of the stellar halo of NGC 5128 based on numerical 
simulations of elliptical galaxy formation. We here adopt the 
basic "merger" approach (Toomre 1977) in which elliptical 
galaxies are proposed to be formed by major merging of two 
spiral galaxies. Within the context of this assumption, we 
investigate the final MDF of the outer halo component of 
the merger remnant (i.e., the giant E galaxy) as well as its 
dependence on the input parameters of the progenitor discs 
(e.g., the MDFs of discs, bulge-to-disc-ratio, and halo mass 
fraction of the disc). Comparing the simulated MDFs with 
the observed ones for NGC 5128, we discuss the following 
points: (1) how the dynamics of galaxy merging is impor- 
tant for the formation of the gE stellar halo, (2) how the 
initial MDF of a disc (or bulge) in a merger controls the final 
MDF of the stellar halo, (3) whether the bulge-to-disc ratio 
of a merger progenitor spiral is important for the determi- 
nation of the stellar halo's MDF, (4) what initial conditions 
of galax;y mergers in the present simulations can best give 
the MDF similar to the observed one, (5) the relevance of 
the MDF for the globular cluster system, which is different 
from the field-halo stars, and (6) whether, on the grounds of 
the merger picture, we should always expect the MDFs of 
stellar haloes in elliptical galaxies (E) to be systematically 
different (more metal-rich) from those in late-type spirals 
(Sp). 

The plan of the paper is as follows: In the next section, 
we describe our numerical models for stellar halo formation 
in galaxy mergers. In §3, we present the numerical results 
mainly on the final MDFs of merger remnants (i.e., elliptical 
galaxies) for variously different merger models. In §4, we pre- 
dict a possible relationship between MDFs of stellar haloes 
and physical properties of their host elliptical galaxies. In 
this section, we also discuss the origin of M3rs metal-rich 
stellar halo. We summarize our conclusions in §5. 



2 THE MODEL 

2.1 Merger model parameters 

In the model calculations, we numerically investigate dy- 
namical evolution of major galaxy mergers between spiral 
galaxies whose primary components arc dark matter haloes, 
stellar discs, stellar haloes, and globular clusters (GCs). In 
the present study, we do not include gaseous components 
(thus no star formation) in the simulations so that all of the 
halo stars in the merger remnant come from the pre-existing 
stellar components of the progenitor spirals. The density 
profiles of dark matter and disc components in a spiral fol- 
low the Fall-Efstathiou (1980) model. The total mass and 
the size of a disc are assumed to be Afd and R^, respectively. 
All masses and lengths are measured in units of Md and Rd 
unless otherwise specified. Velocity and time are measured 
in units oi v = [GMd/RdY''^ and idyn = (i^d/GMd)^/^ re- 
spectively, and the gravitational constant G, is assumed to 
be 1.0. If we adopt Md = 6.0 x 10^" Mq and Rd = 17.5 
kpc £is fiducial values, then v = 1.21 x 10^ km/s and tdyn 
= 1.41 X 10® yr. 

In the disc model, the rotation curve becomes nearly 



© 1994 RAS, MNRAS 000, 1-?? 



Origin of the metallicity distribution of the NGC 5128 stellar halo 3 



flat at i? = 0.35 (where R is distance from the centre of the 
disc) with the niaxinmm rotational velocity Vm = 1-8 in our 
units (220 km s^^). The corresponding total mass of the dark 
matter halo (within 1.5iid) is 4.0 in our units. The veloc- 
ity dispersion of the dark halo component at a given point is 
set to be isotropic and given according to the virial theorem. 
The radial {R) and vertical {Z) density profiles of the disc 
are assumed to be proportional to exp(— i?/iio) with scale 
length Rq = 0.2 and to sech^ {Z / Z{)) with scale length Z{) = 
0.04 in our units, respectively. In addition to the rotational 
velocity caused by the gravitational field of disc, bulge, and 
dark halo components, the initial radial and azimuthal veloc- 
ity dispersions are assigned to the disc component according 
to the epicyclic theory with Toomrc's parameter Q — 1.5. 
The vertical velocity dispersion at given radius is set to be 
0.5 times as large as the radial velocity dispersion at that 
point, as is consistent with the observed trend of the Milky 
Way (e.g., Wielen 1977). 

The density profile of the central bulge with a mass of 
Mb and effective radius of Re (in our units) is represented by 
the de Vaucouleurs R^^* law. We mainly investigate two dif- 
ferent bulge models: (a) a model with Mb = 0.2 (in our units, 
corresponding to 1.2 x 10^" M©) and Rc — 0.04 (in our units 
corresponding to 0.7kpc) and (b) a model with Mb = 0.6 (3.6 
X 10^° Mo) and 0.11 (2kpc). These two cases are intended to 
correspond roughly to the Milky Way and to M31, respec- 
tively. Each spiral is assumed to have an extended stellar 
halo with a density distribution of p(r) ~ r~^'^, where r is 
the distance from the center of the spiral. The mass fraction 
of the stellar halo is represented by /h and considered to be 
a free parameter in the present study. 

We also include a population of coUisionless particles 
to represent a part of the globular cluster system in each 
galaxy. The Milky Way and M31 are observed to have 160 
± 20 and 400 d= 55 globular clusters, respectively (van 
den Bergh 1999). Wc therefore assume each spiral to have 
200 old, metal-poor globular clusters with a number den- 
sity distribution like that of the Milky Way system (i.e., 
p(r)-r-3«). 

In total, the numbers of particles used for modeling the 
merger arc 20000 for the dark matter haloes, 20000 for the 
discs, 12000 (4000) for bulges (The Galactic bulge), 4000 for 
the stellar haloes, and 400 for GCs. 

The orbits of the two spirals in any given run of the 
model are set to be initially in the xy plane by definition, 
and the distance between the centre of mass of the two spi- 
rals (rin) is 6Rd (105 kpc). The pericentre distance (rp) and 
the orbital eccentricity (cp) are assumed to be free parame- 
ters which control the orbital angular momentum and energy 
of the merging galaxy. Wc show here the particular cases for 
Bp =1 and rp = _Rd (17.5 kpc), because the dynamics of halo 
formation do not depend strongly on these merger parame- 
ters. The spin of each spiral in a merging galaxy is specified 
by two angles 9i and (j>i (in units of degrees), where the 
suffix i is used to identify each galaxy. Here, 9i is the angle 
between the z axis and the vector of the angular momentum 
of the disc, and 4'i is the azimuthal angle measured from x 
axis to the projection of the angular momentum vector of 
the disc onto the xy plane. We specifically investigate the 
following three models with different disc inclinations with 
respect to the orbital plane: (a) a prograde-prograde model 
represented by "PP" with Oi =0,02= 30, 0i = 0, and ^2 = 



0, (b) a prograde-retrogradc ("PR") with Oi = 30, 62 = 240, 
(pi = 90, and (f>2 = 0, (c) a retrograde-retrograde ("RR") 
with Oi = 180, 02 = 210, 01 = 0, and 02 = 0. All calcula- 
tions related to the above coUisionless evolution have been 
carried out on the GRAPE board (Sugimoto et al. 1990) at 
the Astronomical Institute of Tohoku University. 



2.2 MDFs of spirals 

A key point of these simulations is to trace out what gener- 
ates the MDF of the merger product halo. We thus assume 
that each stellar component within the progenitor spirals has 
its own characteristic MDF which does not change during 
a simulation of a galaxy merger. The initial MDFs for each 
of the stellar components (i.e., bulge, disc, halo, and GC) in 
a merger progenitor spiral are assumed to be different from 
one another, and are defined to mimic the observed ones in 
the Milky Way. 



2.2.1 Disk 

Recent photometric and spectroscopic observations of the 

Galactic stellar and gaseous components and gaseous com- 
ponents of external galaxies have demonstrated that disc 
galaxies have metallicities that depend ou location within 
the disc, i.e., metallicity gradients (e.g., Zaritsky, Kennicutt, 
and Huchra 1994; Friel 1995). Considering these observa- 
tions, wc allocate metallicity to each disc star according to 
its initial position: at r — R, where r (R) is the projected 
distance (in units of kpc) from the center of the disc, the 
metallicity of the star is given as: 

[m/H]^^j, = [m/H]^ + Od x R. (1) 

We adopt two plausible values for the slope Od- (a) —0.091 
from Friel (1995) in which the Galactic stellar metallicity 
gradient is estimated from Galactic open clusters, and (b) 
—0.056 from Zaritsky ct al. (1994) in which gaseous metallic- 
ity gradients of disc galaxies are estimated from HII regions 
of these. (It should be noted here that we do not have any 
observational results which have directly measured stellar 
metallicity gradients of external disc galaxies.) The models 
are labeled £is "OC"-type (those with the Friel et al. gradi- 
ent defined by open cluster metallicities) and "GAS" -type 
(those with the Zaritsky et al. gradient based on nebular and 
gaseous abundances). The value of the zeropoint [m/H]^ 
is chosen such that the mean metallicity of a disc (repre- 
sented by [m/H]^ = —0.15 for the fiducial model) can be 
different for different models. 

We comment that the steeper "GC" slope might a pri- 
ori be considered the preferable one, since it is based on the 
old- to intermediate-age stars which represent the bulk of 
the stars in the disc, and these stars are likely to be most 
relevant to mergers at the earlier epochs we are attempt- 
ing to model here. By contrast, the "GAS" slope represents 
the youngest stars and present-day gas, which are a much 
smaller fraction of the disc mass. As will be seen below, aa 
turns out to be perhaps the single most important parame- 
ter for understanding the outcome of the post-merger halo. 



© 1994 RAS, MNRAS 000, 1-?? 



4 K. Bekki, William E. Harris, Gretchen L. H. Harris 



2.2.2 Bulge 

The Galactic bulge is also observed to have a metallicity 
gradient (e.g., McWilliam & Rich 1994; Wyse at al. 1997; 
Frogel at al. 1999). Tharafora in the models wa assign tha 
metallicity of a bulge star aX, r = R, where r (R) is the 
projected distance (in units of kpc) of the star from the 
center of the bulge, to bo: 

[m/iil,,=n = [^/iil.r=o + «b X R. (2) 

We adopt the metallicity gradient from Frogel et al. (2000) 
in which the slope Ob = —0.43. For the model with Afi, 
= 0.2, we use [m/H],^ = 0.034. Since there is an ob- 
served relation between luminosity (L) and metallicity (Z) 
for spheroidal galaxies (L oc Z" -*; Mould 1984), we allocate 
a larger metallicity of [m/HJj^ to the model with a larger 
bulge. For example, the value of [m/H]|^ in the model 
with Mb = 0.6 is set to be 0.19. 

2.2.3 Halo and GCs 

There is no clear evidence of radial metallicity gradients in 
the Galactic halo stars or halo globular cluster system (e.g., 
Searle & Zinn 1978; Chiba & Yoshii 1998). We accordingly 
do not allocate a metallicity to each halo star (GC) based 
on its initial position. Instead, we assume that the halo stars 
(GCs) in a model have a Gaussian metallicity distribution 
with a mean of [m/H]j^ ([ni/H]g;.) = —1.4 and dispersion of 
cr([m/H]) = 0.3 (and truncated at 4(j). The adopted value of 
cr([m/H]) = 0.3 has already been demonstrated to be a rear 
sonable value in explaining the observed MDF of the Galac- 
tic stellar halo and GCs (e.g.. Cote et al. 2000). We set the 
value of [m/H]g^ to be —1.6 for all models. 

Most large galaxies also have GCs which are more 
metal-rich and usually more concentrated to the galaxy cen- 
ter than the metal-poor halo clusters. In the Milky Way and 
M31, these more metal-rich clusters (which have a mean 
[m/H] ~ —0.5 and higher overall rotation speed around the 
galaxy center) form about a quarter to a third of the whole 
GC population in both M31 and the Milky Way, and can 
plausibly bo interpreted as a "bulge" component (see Per- 
rett et al. 2002; Harris 2001; Minniti 1996). The MDF for 
the entire globular cluster system in these and many other 
large galaxies can usually be described as having a bimodal 
form (e.g. Larsen et al. 2001), often with a clear separation 
between the halo and bulge components. We do not explic- 
itly include any bulge GCs in the simulations, because (as 
will be seen below) they behave similarly to the bulge stars 
during the merger, and at this stage we are interested pri- 
marily in the formation of the elliptical halo. We will discuss 
the GC results further in §3.2. 

The MDFs for the various components of a progenitor 
spiral are given in Fig. 1 for the fiducial model (described 
in detail later) with the "OC" type metallicity gradient, 
[m/H]^ = -0.15, Mb = 0.2 (The Galactic bulge model), 
and [m/H]i^ = -1.4. 

2.3 Main points of analysis 

For the NGC 5128 halo, MDFs have been obtained at pro- 
jected distances of 8 kpc, 21 kpc, and 31 kpc from the galaxy 
center (HHP, HHOO, HH02). The two outer-halo locations 




I I I I I I I I I I I I I I I L 

-3 -2 -1 

[m/H] 



Figure 1. Upper: The MDFs of each of stellar component, disc 
(thick solid), halo (thin solid), and bulge (dotted) in a spiral 
galaxy for the fiducial model Ml. "OC"-type metallicity gradi- 
ent and a bulge mass (Mb = 0.2) are assumed in this model. 
Lower: The MDFs of the stellar halo (solid) and globular clus- 
ters (represented by GC, dotted) in a spiral galaxy for the fiducial 
model Ml. Here [m/H]^ = —1.4 and [m/H]^^ = —1.6 are assumed 
in this model. 
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Figure 2. Final mass distribution of stellar components that 
are initially disc stars (upper left), halo ones (upper right), GCs 

(lower left), and bulge stars (lower right) in the elliptical galaxy 
formed by major merging in the fiducial model Ml at T = 20 
in our units (2.8 Gyr). One frame measures 77 kpc on a side. 
Note that the halo region (i? > 0.5) of this elliptical galaxy is 
dominated by disc stars. 
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Figure 3. Radial dependence of mass fraction of each stellar 
component, disc (solid), halo (dotted), and bulge stars (dashed) 
in the merger remnant for the fiducial model Ml. Note that more 
than 90 % of stars with R > 0.5 (regarded as "stellar halo" of the 
elliptical in the present study) are initially stars in the disc of the 
merger progenitor spirals. 



(21 and 31 kpc) turned out to give identical results for the 
MDF. For comparison with the models, we investigate the 
model MDFs of two halo regions which correspond closely 
with the these inner and outer regions, one at 0.5 < J? < 
1.0 (referred to as "inner halo") and the other at 1.0 < 7? 
< 2.0 ("outer halo") in each merger model. The simulated 
MDF of the "inner halo" and that of the "outer halo" are 
compared with the observed MDF derived for the 8 kpc field 
and with that derived for the combined 21 and 31 kpc fields. 
Although wc investigated ~ 20 models, we here present only 
the results of the models which are not obviously inconsistent 
with observational results. That is, we attempt to outline the 
range of model parameter space which leads to realistic re- 
sults. For example, our simulations have shown that, if the 
haJo mass fraction of the merger progenitor spirals is greater 
than 0.2, the resultant MDF has far too many metal-poor 
stars to agree with the observed MDF. Equally, if the ini- 
tial disc's mean metallicity in a merger is much larger than 
[m/H] = 0.0, the peak value of the resultant MDF is too 
metal-rich to be consistent with observations. We will not 
show the results of these obviously inconsistent models here 
but, rather, we present the results of the more realistic mod- 
els in which the mean disc metallicity ([m/H]^ ) is equal to 
-0.15 and 0.0, and in which the halo mass fraction is no more 
than 0.2. 

Below, we describe the results of 10 models and in Table 
1 summarize their input parameters: model number (column 
1), bulge mass fraction Mb (column 2), metallicity gradi- 
ent type (3), mean metallicity of disc in [m/H]^ (4), orbital 
type of galaxy merging (5), halo mass fraction /h (6), and 
initial halo mean metallicity [m/H]j^ (7). For convenience 
we refer to the model Ml as the fiducial model. It shows 
typical behavior of stellar halo formation in major merg- 
ing and its MDF is reasonably consistent with the observa- 
tions of the NGC 5128 stellar halo. Accordingly we describe 
mainly the results of this model in the following sections. 
The MDFs of this model for each component (disc, bulge, 
halo, and GCs) are summarized in Table 1. The metallic- 
ity gradient types "OC" and "GAS" in the third column 
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Figure 4. Morphological evolution of disc stars that become halo 
stars with R > I (where R is the projected distance from the 
center of the merger in our units) at T = 20 (2.8 Gyr) in the 
fiducial merger model Ml projected onto the x-y plane. The scale 
is given in our units, and accordingly each frame measures 175 kpc 
on a side. Note that the outer halo stars of the merger remnant 
are initially located in the outer part of the progenitor discs. Note 
also that tidal stripping (or angular momentum redistribution) of 
disc stars during major galaxy merging is a key physical process 
of the stellar halo formation in the elliptical galaxy. The six boxes 
show this evolution at times of T=0,4,8,12,16,and 20. 



represent the two choices for disc metallicity gradient de- 
scribed above. Orbital type labels "PP" , "PR" and "RR" 
represent a nearly prograde-prograde orbital configuration, a 
prograde-retrogradc one, and a nearly retrograde-retrograde 
one. The structural properties of merger remnants are briefly 
summarized in Appendix A. 



3 RESULTS 

3.1 Dynamics of stellar halo formation and MDFs 
of ellipticals 

Figs. 2 and 3 summarize the final mass distributions of the 
stellar components that are initially within discs, haloes, and 
bulges of the merger progenitor spirals for the fiducial mar- 
jor merger model Ml. Since the bulges initially have more 
compact configurations than the discs, only a small fraction 
of bulge stars are tidally stripped to form the inner stellar 
halo with 0.5 < R < 1.0 of the merger remnant, where R is 
the radius from the center of the developed elliptical galaxy 
in units of 7?^- The stars in the discs, on the other hand, 
are very efficiently stripped owing to the strong tidal field 
of major merging (or by angular momentum redistribution 
during merging) so that they end up in both the inner and 
outer halo (1.0 < R < 2.0) of the elliptical galaxy. Although 
the initial stellar haloes of the spirals are also stripped (or 
'puffed up') efficiently to form the inner and the outer stellar 
haloes, the stellar halo of the resulting elliptical is dominated 
by the stars that were initially within the progenitor discs; 
about 95% of the elliptical halo comes from the initial disc 
stars at 1 < i? (For comparison, the outer parts of the spi- 
rals with 1 < i? < 2 are by definition dominated only by 
their metal-poor halo stars). These results imply that major 
mergers can efficiently populate the stellar haloes of their 
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Figure 5. Radial distribution of the surface (number) density of 
stars located in the halo region of the merger model Ml at T = 20 
(solid) and those in the halo region of the progenitor spiral at T = 
(dotted). Here the results are given for the "outer" halo region 
(1 < R < 2). The scale is given in our units, and accordingly the 
outer halo region corresponds to 17.5 < R < 35.0 kpc. 
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product elliptical galaxies by disc stars from within their 
progenitor spirals. 

Wc note also from Fig. 4 that the (outer) stellar halo of 
the elliptical galaxy comes from the outer part of the merger 
progenitor disc, because the tidal stripping (or angular mo- 
mentum redistribution) is more efficient for the outer part of 
the disc than for the inner part. Furthermore, Fig. 4 shows 
that, although the violent dynamical relaxation of major 
merging results in strong mixing of stars of the two discs, 
we can clearly see the inhomogeneous distribution of halo 
stars at 2 < J? (e.g., the upper right part of the halo ap- 
pears to be denser than the lower left). This implies that 
if we can observe the global stellar halo or surface intensity 
distribution of an elliptical galaxy, it may show a certain 
degree of inhomogeneity. 

As is shown in Fig. 5, the surface density of the merger 
halo is greater by roughly 2 orders of magnitude at a given 
radius than the initial spiral halo owing to the stripped 
disc stars. This result suggests that the surface brightness 
of stellar haloes in elliptical galaxies will be systematically 
higher than for haloes in spiral galaxies if most ellipticals 
are formed by major merging. However, the slope of the 
halo light distribution does not change significantly com- 
pared with that of the initial spiral halo. 

A significant result of the halo formation by stripping 
of disc stars is that the MDF of the stellar halo of the ellip- 
tical galaxy becomes remarkably different from that of the 
progenitor spirals. We show this, and the comparison with 
the data, in Figures 6 and 7. By comparing Fig.l and Fig. 6, 
we also note that the difference in MDFs between the halo 
and the disc (spheroidal) components in this model becomes 
much smaller after merging. This implies that the MDFs of 
the stellar haloes in elliptical galaxies are less different from 
those of the main spheroidal bodies than they are in spirals. 

As is shown in the lower panel of Fig. 6, the stellar halo 
of the emergent E galaxy is a genuine composite popula- 
tion in this formation scheme. The relatively metal-rich part 
which is the vast majority ([m/H] > —1) come from the pro- 
genitor discs, whereas the thin metal-poor part comes from 
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Figure 6. Upper: MDFs derived for observations (solid) and the 
simulation of model Ml (dotted). The MDF for observations is 
constructed from halo stars with the projected distance of 21 and 
31 kpc in NGC 5128 halo fields. The MDF for the simulation is 
constructed from the simulated stars with 1 < i? < 2 (17.5 kpc < 
R < 35 kpc), where R is the distance of a star from the center of 
the merger remnant. Note that the simulated MDF shows a peak 
around [m/H] ~ —0.4, which is similar to the observed peak. 
Lower: MDFs for stars that are initially disc stars (solid) and 
halo ones (dotted) in the simulation of model Ml. These MDFs 
are also estimated for the outer halo region with 1 < ij < 2. Bulge 
stars do not contribute to the MDF in this outer halo at all. 



the progenitor haloes. A very small fraction at the metal- rich 
end is from the original bulges. 

Fig. 7, showing the inner halo (0.5 < R < 1.0) of the 
elliptical galaxy, has an MDF which is shifted even further 
to the metal-rich end, with a broad peak at [m/H] ~ —0.4 
and an extension to well above Solar metallicity, and an 
even sparser low-metallicity tail than the outer halo. Here, 
a higher fraction of the stars are from further in within the 
progenitor discs (~ 97 %), as well as from the initial bulges 
(~ 2 %). The simulated MDF is not as broad as the ob- 
served inner halo MDF and, in particular, does not extend 
as far into the metal-rich regime. This may be significant, 
particularly since the observed data are likely to be more 
incomplete for the reddest stars (HH02). Thus the principal 
driver of the observed metallicity gradient in the resulting E 
galaxy halo is the presence of the initial metallicity gradients 
in the disc. Fig. 8 demonstrates that the MDF in the very 
outer part of the halo {R > 2, corresponding to i? > 35 kpc) 
shows a much wider peak around — 1 < [m/H] < —0.4 and a 
larger fraction of metal-poor ([m/H] < —1.0) stars compared 
with the MDFs for the inner and the outer haloes. This is 
mainly because the very outer halo with R > 2 consists of 
more metal-poor stars coming from the extreme outer parts 
of the merging discs. The predicted difference in the MDFs 
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Figure 7. The same as the Figure 6 but for tlie inner halo region 
with 0.5 < R < 1.0. The MDF for observations is constructed 
from halo stars with the projected distance of 8 kpc in NGC 5128 
halo fields. In the lower panel, MDFs for stars that are initially 
disc stars (solid) and halo ones (dotted) , and bulge ones (dashed) 
are plotted. 



between the outermost halo (2 < R) and the already ob- 
served inner and outer haloes (i.e., 0.5 < R < 1.0 and 1.0 < 
R < 2.0) can be tested by future observations on the MDF 
of the NGC 5128 stellar halo at projected distances of larger 
than 30 kpc. 

In Fig. 9, we show the dependence of mean metallicity 
on galactocentric distance for the merged elliptical and the 
initial spirals. The outer stellar halo of the model elliptical 
(1 < < 2) does not show any important metallicity gra- 
dient, though its progenitor discs have negative metallicity 
gradients. This is mainly because violent dynamical relax- 
ation during galaxy merging can cause efficient mixing of 
stellar populations (with different metallicity) initially lo- 
cated in different regions of the discs, and this almost com- 
pletely smoothes out the initial disc's metallicity gradient. 
(It should be noted here that if we observed the gradient 
throughout the elliptical, we can still see the negative gradi- 
ent, which becomes steeper farther in to the center.) Fig. 9 
also demonstrates, as indicated above, that the mean metal- 
licity of the stellar halo of the elliptical is much higher (~ a 
factor of 6) than that in the progenitor spiral at all points in 
the outer halo regions. This is remarkably consistent with 
the HST data for the NGC 5128 fields at 21 kpc and 31 
kpc whose MDFs are virtually indistinguishable (HHOO). We 
suggest that the near-zero metallicity gradient of the outer 
stellar halo is an observational test which can assess the va- 
lidity of the merger scenario of elliptical galajcy formation. 

In Fig. 10 we show the distribution of the halo stars on 
the l/y-[m/H] diagram (radial velocity vs. metallicity) for the 
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Figure 8. The same as the Figure 6 but for the very outer halo 
regions with 2.0 < R. The MDF for observations is constructed 
from halo stars with the projected distance of 21 and 31 kpc in 
NGC 5128 halo fields. 
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Figure 9. Metallicity gradient of the stellar halo with 1 < i? < 
2) in the initial spiral (solid) and the merger remnant (dotted) for 
the model Ml. Note that values of log (R/Rd) = 0.08 and 0.25 
correspond to the approximate locations of the observed 21kpc 
and 31 kpc fields. 



progenitor spirals and the elliptical. Only the elliptical halo 
contains stars with both high projected (radial) velocity {Vy 



2 in our units, corresponding to ^ 280 km 



and high 



metallicity ([m/H] > —0.6). Some fraction of these outer 
halo stars with rather high projected velocity also have large 
tangential velocities and thus high angular momentum with 
respect to the elliptical. The origin of these halo stars with 
high angular momentum is closely associated with angular 
momentum transfer, from inside to outside, during major 
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Figure 10. Distribution of halo stars with 1 < R < 2 (outer halo) 
on Vy-[m/H] diagram for the merger progenitor spiral (upper) and 
the merger remnant (lower) in the model Ml. Vy here represents 
the velocity of a star projected onto y-axis of the simulation. Note 
that there are a larger number of halo stars with [m/H] > —0.6 
and the absolute magnitude of Vy larger than 1.8 (250 km s~^) 
only in the merger remnant. A physical interpretation of this is 
given in the text. 
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Figure 11. MDFs of halo stars (solid) and halo GCs (dotted) in 
the merger remnant of the model Ml at T = 20. The MDFs are 
constructed from the stellar components with 0.5 < R (i.e., both 
the inner and the outer stellar components). Note that there is 
a clear difference in MDFs between the halo and the GC. This 
clear different is not seen in the initial spiral's MDF (See Fig. 1) , 
which suggests that dynamical processes of galaxy merging cause 
such a remarkable difference in the MDFs between the halo and 
the GCs in the model Ml. 



in NGC 5128 are beyond current observational capabilities, 
the results in Fig. 10 suggest that, future such observations 
may find a larger fraction of relatively metal-rich outer-halo 
stars with high angular momentum in NGC 5128. 



3.2 Globular clusters 

By comparing Fig. 1 with Fig. 11, we can clearly see that 
the diflerence between the GC and field-star MDFs in this 
merger model becomes much more conspicuous after ellip- 
tical galaxy formation. This is essentially because the GCs 
start out as part of a "hot" , pressure-supported component 
of the spiral haloes and so behave differently from the disc 
stars during the merger. The MDF of the halo clusters thus 
does not change markedly during merging, whereas that of 
the halo can be dramatically changed owing to the addition 
of large numbers of stripped disc stars. In the product ellipti- 
cal, we therefore end up with metal-poor clusters embedded 
in a moderately rnctal-rich and massive halo. 

This dynamically constructed difference in MDFs would 
not be changed if we were to include gas and merger-induced 
formation of metal-rich GCs, since numerical simulations 
have already demonstrated that these would mostly be lo- 
cated in the central region of the merger remnant and there- 
fore would not contribute to the MDF of GCs in the outer 
halo region (Bekki et al. 2002). In addition, the "bulge" GCs 
which were already present in the progenitor spirals (such as 
are present in both the Milky Way and M31) are more con- 
centrated to their parent galaxy centers and would stay in 
the inner regions along with the bulge stars. The significant 
difference in MDFs between the halo stars and the GCs that 
we note in the Ml simulation has been actually observed 
for M31 (Durrell et al. 2001) and NGC 5128 (HHP, HHOl, 
HH02), both of which have large central spheroidal compo- 
nents. A mean offset in integrated color (and by hypothesis 
metallicity) between GCs and field stars is also classically 
observed in many giant ellipticals (e.g. Harris 1991, Brodie 
& Huchra 1991). It has frequently been suggested that the 
redder, more metal-rich clusters arc similar to the bulk of the 
field halo and bulge stars, while the bluer, more metal-poor 
ones have earlier and different origins, such as within spiral 
haloes or dwarf galaxies (e.g. Larsen et al. 2001; Forbes, 
Brodie, & Grillmair 1997; Geisler, Lee, & Kim 1996; Cote, 
Marzke, & West 1998). The Milky Way, an Sbc-type spiral 
with only a small bulge, is not observed to have such a differ- 
ence between halo stars and clusters. We therefore suggest 
that the observed difference in MDFs between stellar haloes 
and GCs in elliptical galaxies or early-type spirals with big 
bulges is consistent with the view that their large spheroidal 
components were formed by past major merger events. We 
discuss this point further in §4. 



galaxy merging; orbital angular momentum of the merging 
spirals is converted into the intrinsic angular momentum of 
such outer halo stars. Based on kinematics and spatial dis- 
tribution of the planetary nebula (PN) system in NGC 5128, 
Hui et al. (1995) found a possible evidence of the PN's rota- 
tion increasing with radius and discussed the origin of this 
elliptical galEixy. Additional new velocity analyses are given 
by Peng et al. (2002). Although the necessary spectroscopic 
data to obtain kinematical properties of individual halo stars 



3.3 Parameter dependences 

Although the models Ml - MIO resemble each other in their 
dynamics of stellar halo formation and the peak value of the 
MDFs of the resulting stellar haloes, the shapes of the MDFs 
depend on the progenitor disc metallicity gradient, bulgc-to- 
disc-ratio, and the mass fraction and the MDFs of the initial 
haloes. In Fig. 12, we illustrate these effects: 

(i) The MDFs show a larger fraction of metal-rich (—0.4 
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Figure 12. Dependence of the simulated MDFs on the initial 
mass fraction of bulge (top), the MDFs of the merger progen- 
itor discs (second from the top), orbital configurations of ma- 
jor merging (second from the bottom), and MDFs of the stellar 
haloes of the merger progenitor spirals (bottom). For compar- 
ison, the MDF observed for the stellar halo regions with pro- 
jected distances of 21 and 31 kpc in NGC 5128 (HHP, HHOO, 
and HH02) is given by a thick solid line in each panel. The re- 
sults are given for (a) the model with Mj, = (solid, M3), (b) 
Mb = 0.2 (dotted, Ml), (c) Mb = 0.6 (dashed, M2), (d) "OC- 
type metallicity gradient (solid. Ml), (e) "Gas"-type one (dotted, 
M4), (f) more metal-rich disc with [m/H]j = (dashed, M5), (g) 
"PR" (prograde-retrogradc) orbital configuration (solid, M6), (h) 
"PP" (prograde-prograde) orbital configuration (dotted, M7), (i) 
"RR" (retrograde-retrograde) orbital configuration (dashed, M8), 
(j) the fiducial set of parameters (solid. Ml), (k) more massive 
(/h = 0.2) stellar halo (dotted, M9), and (1) more metal-poor 
([m/H]b = -2) stellar halo (dashed, MIO). 



< [m/H] < 0) halo stars in the model with the larger bulge- 
to-disc-ratio (M2, Mi, = 0.6). This is because more of the 
metal-rich bulge stars can be more efficiently stripped into 
the stellar halo owing to its larger extension of the bulge in 
the model with the larger bulgc-to-disc-ratio. The MDF do- 
rived in the model with the smaller bulge-to-disc-ratio (M3, 
Mi = 0) seems to be more similar to the observed MDF of 
the NGC 5128 stellar halo, although we stress (see HH02) 
that the data may be somewhat incomplete at the highest 
mctallicitios. Nevertheless, the current comparisons would 
allow NGC 5128 to be formed by major merging between 
later-type spirals. 

(ii) The stellar halo in the model (M4, GAS-typc) which 
was generated from progenitor spirals with the shallower disc 
metallicity gradients (gaseous nebulae abundances; Zaritsky 
et al. 1994) can not match the observations. Such models 
predict far too few stars in the intermediate metallicity range 
( — 1.4 < [m/H] < —0.6). The input discs with steeper metal- 
licity gradients (OC-type) produce a much closer match 
to the observations. The opposite problem is presented by 
model M5, which has a larger initial mean disc metallicity; 
this model ends up predicting far too many halo stars for 
[m/H] > —0.2. These results are the most prominent ones 
of our investigation, and indicate that the detailed shape of 
the observed MDF in NGC 5128 can provide information on 
the mean metallicity and disc gradients in the merger pro- 
genitor spirals (and thus, possibly, on their luminosity due 
to the observed color- magnitude relation for disc galaxies). 

(iii) The shapes of the MDFs of stellar haloes of el- 
liptical galax;ies formed by major merging do not depend 
strongly on the orbital configurations of the merging (M6, 
7, and 8). This probably reflects the fact that the physical 
processes of violent relation and angular momentum redis- 
tribution do not so strongly depend on the orbital parame- 
ters of major merging. The PP model shows a slightly higher 
peak value of the MDF, because a larger amount of disc stars 
can be stripped owing to more efBcient angular momentum 
redistribution process in this model. 

(iv) Lastly, the stellar haloes in the progenitor spirals 
exert noticeable effects on the stellar halo of the merger rem- 
nant at the low-metallicity end ([m/H] < —1.2). The model 
M9 with the larger initial mass fraction of the spiral's halo 
(20 %) produces too many halo stars at —1.8 < [m/H] < 
— 1 and is not consistent with the observations. The model 
haloes with mass fractions of 5% yield much closer matches 
to the data. In short, the model comparisons appear to rule 
out progenitor spirals with very massive, metal-poor haloes. 



4 A NEW CONSTRAINT ON ELLIPTICAL 
GALAXY FORMATION 

Although the observed tight color-magnitude relation of el- 
liptical galaxies has been considered to suggest that elliptical 
galaxies are old, coeval and homogeneous systems passively 
evolving after the single initial burst of star formation asso- 
ciated with dissipative galaxy formation (e.g., Bower, Lucey, 
& Ellis 1992), a growing number of recent observational re- 
sults suggest that elliptical galaxies exhibit a great variety 
of star formation history (e.g., Worthey, Trager, & Faber 
1996). The observed narrowness of the fundamental plane 
of elliptical galaxies implies a rather smaller range of per- 
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mitted dynamical state of the galaxies (Djorgovski & Davis 
1987; Dressier et al. 1987). However, the morphological di- 
chotomy between boxy-disky elliptical galax:ies (Kormendy 
& Bender 1996) and departures of the projected density pro- 
files from the de Vaucoulcurs R^'''^ law (Caon, Capaccioli, 
& D'Onofrio 1993) show a great variety of major stellar or- 
bit families within the galaxies. Elliptical galaxies thus show 
both diversity and uniformity in their chemical, photometric 
and dynamical properties which impose valuable constraints 
on any theoretical models of elliptical galaxy formation. 

The present study has demonstrated that, irrespective 
of bulge-to-disc-ratio and MDF of of progenitor spirals, el- 
liptical galaxies formed by major merging show the ten- 
dency for stellar haloes to be populated predominantly by 
moderately metal-rich stars in the range — 1 < [m/H] < 
which were originally in the spiral discs. Classical mono- 
lithic collapse scenarios (e.g., Larson 1976; Carlberg 1984) 
which have shown strong negative metallicity gradients (i.e., 
more metal-poor in the outer part of a giant elliptical) are 
less likely to yield the moderately metal-rich populations in 
the outer halo regions that we see in our test case of NGC 
5128. Therefore, if future systematic observations on MDFs 
of elliptical galaxies confirm the above generic trend in their 
halo MDFs, it would provide evidence more consistent with 
major merging as a promising scenario of elliptical galaxy 
formation. 

Spiral galaxies exhibit a color-magnitude (CM) relation 
with a slope steeper than those of elliptical galaxies (de 
Grijs & Peletier 1999). If the CM relation is due mostly 
to metallicity vs. luminosity (i.e., no significant age effects 
on integrated colors) , then more luminous spirals should be 
more metal-rich. The present numerical studies have demon- 
strated that the mean halo metallicity of an elliptical de- 
pends strongly on the mean metallicity of the merger pro- 
genitor discs, and furthermore that the mass fraction of stel- 
lar haloes is not different between different merger models. 
These observational and theoretical results lead to the pre- 
diction that if the majority of E galaxies were formed by mar 
jor mergers (that is, a single merger event between roughly 
equal progenitors), there should be a strong correlation be- 
tween the color (mean metallicity) of the stellar halo and the 
absolute magnitude of its host elliptical galaxy: the more lu- 
minous (redder) elliptical should show redder halo colors. In 
this picture, if there are very luminous ellipticals with rather 
metal-poor haloes that violate this correlation, they would 
have been formed by a longer series of mergers from smaller, 
more metal-poor progenitors or perhaps by accretion as in 
Cote et al. (1998, 2002). 

4.1 Origin of the M31 stellar halo 

Since Mould & Kristian (1986) first resolved clearly the 
brightest RGB and AGB stars in the M31 halo, many 
observational studies of its MDF have been done (e.g., 
Pritcliet & van den Bergli 1988; Christian & Heasley 1991; 
Davidge 1993; Durrell, Harris, & Pritchet 1994; Couture et 
al. 1995; Holland, Fahlman, & Richer 1996; Rich, Mighell & 
Neil 1996; Reitzel, Guhathakurta, & Gould 1998; Reitzel & 
Cuhathakurta 2001; Durrell, Harris, & Pritchet 2001). One 
of the most remarkable results of these studies is that the 
derived MDF of the M31 halo is very different from that of 
the Milky Way halo. All the studies listed above find that 



the M31 halo is dominated by moderately high-metallicity 
([m/H] ~ —0.5) populations at a very similar level to what 
we find in NGC 5128. Durrell, Harris, & Pritchet (2001) 
point out that the bulk of the M31 halo stars are much 
less likely to come from dwarf galaxy populations, just as 
for NGC 5128. Cote et al. (2000) demonstrated that if the 
M31 halo was formed by dissipationless merging of numer- 
ous subgalactic clumps with a luminosity (L) function of 
slope dn/dL oc L~^'^, the observed MDF could be repro- 
duced. However, it is not clear a priori why the MDFs of the 
Milky Way and M31 should be so different in the same Local 
Group environment with the plausible assumption that the 
surrounding pregalactic population of clumps would be sim- 
ilar for both. Ibata et al. (2001) suggest that the bulk of the 
M31 halo may have arisen from a single event, namely can- 
nibalization of stars from the similarly metal-rich satellite 
M32. 

We here propose that the origin of the M31 stellar halo 
can be closely associated with M31's bulge formation via 
major merging of two discs at an early epoch (see also Free- 
man 1999). In this scenario, the metal-rich halo stars were 
initially within the outer parts of the merger progenitor spi- 
rals and tidally stripped during the M31's bulge formation. 
This merger scenario can give a plausible explanation as to 
why the MDFs in the haloes of M31 and NGC 5128 re- 
semble each other so strongly (Harris & Harris 2001). If 
the observed difference in the metallicity of the dominant 
metal-rich population between these two can reflect the dif- 
ference in mean metallicity of merger progenitor spirals be- 
tween these two, we can estimate the luminosity difference 
between the progenitor spiral of M3rs bulge and that of 
NGC 5128. The metallicity of the dominant halo popula- 
tion is [m/H] ~ -0.5 for M31 (Durrell, Harris, & Pritchet 
2001). and -0.4 for NGC 5128 (HHP; HHOO, 02). By us- 
ing the observed luminosity-metallicity relation Z oc L"'^ 
(Mould 1984), we can expect that the possible luminosity 
difference of merger progenitor spirals between these two 
big spheroidal components is a factor of ~ 2. This difference 
is appreciably smaller than the observed luminosity differ- 
ence (a factor of ~ 4) between M31's bulge with Ly ~ 1.1 
X 10^° 1/0 corresponding to about 40 % of total stellar light 
of M31 (van den Bergh 2000) and NGC 5128 with Ly ^ 4.7 
X 10^° Lq (Hui et al. 1995). 

Freeman (1999) suggested that the big M31 bulge could 
be formed by a past merger event on the basis of M31's 
halo structural properties. He also pointed out that if M31's 
bulge was formed by merging, its outer stellar halo com- 
ponents should show significant rotation owing to efficient 
angular momentum transfer during merging. Our present 
simulations have confirmed that the outer stellar halo com- 
ponents clearly show rotation after major merging. We note 
that Perrett et al. (2002) have indeed found a strong rotation 
signature for both types of globular clusters in the M31 halo 
(160 km s~^ for the metal-rich population, 130 km s~^ for 
the metal-poor population) , extending outward to ~ 15 — 20 
kpc. 

Big spheroidal bulges (like M31's) are suggested to be 
physically different from box- and peanut-shaped bulges 
in spirals (e.g., Wyse et al 1997 and references therein): 
Spheroidal bulges are "true" bulges similar to giant el- 
lipticals whereas the box- and peanut-shaped bulges are 
edge-on bars formed through dynamical instability of discs 
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(e.g., Combes & Sanders 1981). We suggest that if "true" 
spheroidal bulges are formed by merging (rather than disc 
instability), spirals with such big bulges are more likely to 
show halo MDFs similar to that of M31 rather than to that 
of the Milky Way. It would be an interesting future study 
to correlate the MDFs of stellar haloes with bulge shapes in 
different Hubble types of galaxies. 

4.2 Discussion and Remaining Problems 

We close by briefly mentioning alternate models of halo for- 
mation applicable to NGC 5128 and other ellipticals, along 
with difficulties still faced by most interpretations. 

E galaxies may also form by longer series of mergers (ac- 
cretion) of smaller subgalactic clumps. The dissipationless 
accretion of dwarf galaxies with appropriate chosen lumi- 
nosity functions can be also responsible for MDFs with rel- 
atively metal- rich peaks (e.g., Cote ct al 1998, 2000, 2002). 
As also noted earlier, Beasley et al. (2002a,b) take the alter- 
nate route of hierarchical merging of pregalactic gas clouds 
to model the NGC 5128 MDF, with first-order matches pos- 
sible in a wide variety of circumstances. Without repeating 
the extensive analysis in these papers, we note only that 
their work demonstrates that a unique physical interpreta- 
tion for the observed MDFs may still not be in hand. Even- 
tually, all the observational parameters including kinematics 
and halo/bulge structures as well as MDFs must be brought 
to bear on the question. 

There are, however, two interesting interpretive prob- 
lems which the globular cluster populations confront us 
with. First is the sheer size of the cluster population in 
NGC 5128. Coincidentally, the luminosity of NGC 5128 
(My ~ —22.0) is similar to what we would get if we sim- 
ply combined two galaxies like M31, or M31 plus the Milky 
Way (Pritchet & van den Bergh 1999). Thus a simple merger 
of the two would give an elliptical of about the same size 
as NGC 5128 and with - 600 - 800 clusters. A GC popu- 
lation of that size is two to three times smaller than the 
~ 1500 - 2000 globular clusters that are in NGC 5128 
(G.Harris et al. 1984). That is, NGC 5128 has a GC specific 
frequency significantly higher than normal spiral galaxies. 
Contrary to many comments in the literature, this classic 
"specific frequency problem" cannot be solved by invoking 
the presence of gas and then forming large numbers of clus- 
ters during the merger. This is because field stars will also 
form from the gas, and the ratio of field stars to GCs will stay 
at much the same level. See the discussion of Harris (2001), 
who notes that to produce noticeable changes in the spe- 
cific frequency, it is necessary to have both extremely large 
amounts of gas and extremely high cluster formation effi- 
ciency. Under such conditions, we would no longer be talking 
about a simple disc/disc merger but rather an earlier epoch 
of primary galaxy formation from hierarchical merging of 
gas clouds. 

A second concern relates to the spatial distributions of 
the GCs in the metal-poor and metal-rich subgroups that 
many or most giant ellipticals contain. The GCs in the 
metal-rich mode are not found only in the bulge region as 
noted above, but also extend outward far into the halo in sig- 
nificant numbers (e.g., G.Harris et al. 1992; Lee & Geisler 
1993; Geisler, Lee, & Kim 1997). Such metal- rich clusters 
are also present in the outer halo of M31 (Perrett et al. 



2002). Within the context of the merger scenario we have 
tested here, it is not clear how these metal-rich, outer-halo 
clusters arise. Clusters formed from appropriately metal-rich 
gas during the merger are likely to build up in the central 
regions and add to the bulge population. Metal-rich globular 
clusters that are already distributed widely throughout the 
discs or haloes of the progenitor spirals would end up pop- 
ulating the halo of the emergent elliptical. Such clusters are 
present to some extent in M31, but not in the Milky Way. In 
summary, it is still difficult for any one model to correctly 
predict the total ensemble of features of the giant ellipticals, 
which may have formed by a combination of routes. 



5 CONCLUSIONS 

We have numerically investigated MDFs of stellar haloes of 
elliptical galaxies formed by major merging in order to elu- 
cidate the origin of the observed MDF of NGC 5128. We 
mainly investigated the MDF of halo stars with R > 
(where R and Rd are the radius from the center of a merger 
remnant and the initial size of the merger progenitor spi- 
ral) in the merger remnant for several models with different 
physical parameters of the progenitor spirals. Our numerical 
study clearly demonstrates the importance of dynamical pro- 
cesses of galaxy merging (i.e., tidal stripping of metal-poor 
disc stars) in determining MDFs of stellar haloes of elliptical 
galaxies. We summarize our principal results as follows. 

(1) The stellar halo of an elliptical galaxy formed by 
major merging is predominantly populated by moderately 
metal-rich stars with [m/H] ^ —0.4, which come from the 
outer parts of the discs of the merging spirals. Furthermore, 
the fraction of metal-poor stars with [m/H] < —1.0 in the 
halo is very small 17 %). The MDF of the merger remnant 
does not show a strong radial dependence for 0.5 < R < 
2.0 (in our units). These results arc broadly consistent with 
recent -ffST observations on MDFs of the stellar halo of NGC 
5128. 

(2) Irrespective of the physical parameters of merger 
progenitor discs, the MDFs of stellar haloes of merger rem- 
nants show both the dominant moderately metal-rich popu- 
lations and the minor metal-poor ones. Thus future system- 
atic observations on MDFs of elliptical galaxies can assess 
the relative importance of major galaxy merging in the for- 
mation of elliptical galaxies by confirming the above generic 
trend of MDFs. 

(3) The MDFs of elliptical galaxies formed by major 
merging depend strongly on the initial MDFs of merger pro- 
genitor discs, because halo stars are dominated by those 
which are initially in the outer part of the discs and tidally 
stripped during merging. 

(4) MDFs of stellar haloes do not depend so strongly on 
the MDFs of the progenitor spiral bulges, essentially because 
bulge stars are not be tidally stripped so efficiently during 
merging owing to bulge's compact configuration. However, 
the shapes of MDFs of a merger remnant at —0.2 < [m/H] 
< can be affected by the bulge-to-disc-ratio of the merger 
progenitor spiral. 

(5) The shape of the halo MDF of a merger remnant at 
[m/H] < —1.0 depends on the mass fraction and the mean 
metallicity of the progenitor spiral's halo. Therefore the ob- 
served shape of the halo MDF at [m/H] < -1.0 in NGC 
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5128 can give some constraints on the stellar populations of 
the progenitor spiral's stellar halo. 

(6) The observed similarity in MDFs of stellar haloes 
between M31 and NGC 5128 can be naturally explained, if 
the M31 bulge was also formed by major merging between 
two spirals with the masses (a factor of 2 — 4) smaller than 
those that are assumed to merge to from NGC 5128 in the 
present study. The observed markable difference in MDFs 
between stellar haloes and GCs for NGC 5128 and M31 can 
be also explained by the merger scenario. 

(7) Our simulations predict that if most elliptical galax- 
ies are formed by major merging, the photometric properties 
(e.g., colors and magnitudes) of a stellar halo in an elliptical 
galaxy can correlate with those of the galaxy (i.e., redder 
Es have redder stellar haloes). Also our study predicts that 
there can be CM relations in stellar haloes of elliptical galax- 
ies. 
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Figure 13. Surface density profile of the model Ml. The profile 
of the simulated elliptical galaxy is well described as the B}/'^ 
law (de Vaucouleurs law) with the effective radius (He) of 0.26 
{Re^^^ = 0.71) in our units (corresponding to 4.6 kpc) in this 
model. This simulated profile is consistent with the observed one 
by van den Bergh (1976) for NGC 5128. 



Vaucouleurs's law. The effective radius of the simulated el- 
liptical galaxy in each model is 4.6 (Ml), 4.6 (M2), 5.3 (M3), 
6.1 (M7), and 7.0 kpc (M8). These derived values are not 
significantly different from the observed one. As is shown in 
Fig. 13, the simulated radial surface density profile in the 
model Ml follows de Vaucouleurs's law. The radial surface 
density profiles in other models can be also described as de 
Vaucouleurs's ones, which is consistent with the above obser- 
vation (2). These comparison thus suggest that our model 
can reproduce the observed structural properties of NGC 
5128 reasonably well. 
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APPENDIX A: LIGHT PROFILES OF MERGER 
REMNANTS 

Although the main purpose of this paper is not to discuss 
structural properties of NGC 5128, it is important for the 
present study to confirm whether our numerical simulations 
can reproduce not only the MDF observed in the NGC 5128 
stellar halo but also its global structural properties. In the 
following, we assume that (1) the distance (d) of NGC 5128 
is 3.9 Mpc (which is the same as that adopted in HHOO 
and HH02) and (2) the simulated density profile corresponds 
to the light profile (i.e., mass-to-light-ratio is constant for 
the entire region of a simulated elliptical galaxy). We then 
compare the observed radial brightness profile of NGC 5128 
by van den Bergh (1976) with the simulated ones in models 
with initial dynamical conditions different with one another 
(i.e., Ml, 2, 3, 7, and 8). 

van den Bergh (1976) found that (1) the effective ra- 
dius of NGC 5128 is about 5!5 or ~ 6.2 kpc for d = 3.9 
Mpc and (2) the radial brightness distribution follows de 
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Table Al. Model parameters 
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